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Abstract 

I adapted a model recently introduced in the context of seismic phenom- 
ena, to study creep rupture of materials. It consists of linear elastic fibers 
that interact in an equal load sharing scheme, complemented with a local 
viscoelastic relaxation mechanism. The model correctly describes the three 
stages of the creep process, namely an initial Andrade regime of creep relax- 
ation, an intermediate regime of rather constant creep rate, and a tertiary 
regime of accelerated creep towards final failure of the sample. In the tertiary 
regime creep rate follows the experimentally observed one over time-to-failure 
dependence. The time of minimum strain rate is systematically observed to 
be about 60-65 % of the time to failure, in accordance with experimental 
observations. In addition, burst size statistics of breaking events display a 
—3/2 power law for events close to the time of failure, and a steeper decay 
for the all-time distribution. Statistics of interevent times shows a tendency 
of the events to cluster temporarily. This behavior should be observable in 
acoustic emission experiments. 

Keywords: Creep, Stress relaxation, Viscoelastic fiber bundle models, 
Time-to-failure singularity, Numerical algorithms 



1. Introduction 

Creep is the time deformation of a material under the application of a 
mechanical load that is below its overall mechanical strength. In many cases, 
creep can lead to the failure of the material after a finite time after the ap- 
plication of the load, and this kind of creep rupture can be one of the mayor 
concerns in mechanical design. Creep rupture follows a sequence of three 
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temporal stages before complete failure occurs. Immediately after the load is 
applied, there is a so-called primary creep regime in which the strain rate in 
the mat erial decays as a function of time. This is known as the Andrade creep 
regime (Undradef Il91(ih . The second stage is one of almost constant strain 
rate, that gives way eventually to the tertiary regime of accelerated creep 
that leads to complete material collapse. In the Andrade regime and the ter- 
tiary creep regime the strain rate is observed to follow power laws in time. In 
the Andrade regime, strain rate e decreases as a negative power of the time 
elapsed since load application e ~ t~ p , whereas in the tertiary creep regime, 
strain rate increases as a negative power of the time remaining to failure, i.e. 
i ~ (tf — t)~ p , where tf is the time at which failure is observed. Experimen- 
tal values typi c ally ob served for the exponents p and p' are both close to 1 



( INechad at al.l . l2005al Jbl). The time to failure tf depends strongly on the ap- 
plied load a. Actually, if a is below some threshold value a c , failure does not 
occur, and only the initial Andrade regime towards a finite deformation state 
is observed. Time to failure is expected to diverge continuously as a — > <r+, 
however, the precise experimental determination of this dependence suffers of 
the problem that the determination of tf is typically destructive, and differ- 
ent samples have to be used for different values of applied stress, introducing 
sample-to-sample fluctuations that complicate the determination of tf. The 
typical dependence of strain rate vs time displays another interesting fea- 
ture: the time at which strain rate is minimum t m , is typically a conserved 
fraction of the time to fail ure, independe n tly of the level of applied stress, 
and details of the samples. INechad at al.l (l2005al ) find a ratio tf/t m ~ 1.58 
for a set of fifteen samples with different characteristics. In addition to the 
average evolution of the creep process, characterized by the temporal evolu- 
tion of the strain, or the strain rate, there is relevant information that can be 
extracted by analyzing the details of the temporal evolution. In fact, creep 
proceeds typically through bursts of activity in which the strain increases in 
a finite amount, followed by periods of quiescence. This pattern of activity, 
which is reminiscent for instance of the Barkhaussen noise in ferromagnets, 
or earthquakes in seismic processes, can be typically foll owed by recording 
the pattern s of acoustic emission generate d by the sample ( jMaes et al.l . Il998 ; 



Petri at al.l 11994 ; lYagamuchi at al.l . Il99lf ) 



There is a variety of models that have been used to explain part of the 
phenomenology previously described about the creep rupture of materials. 
Most of them are based in the idea of "fiber bundles" , in which the material 
is considered as a set of fibers that fail individually when some threshold 
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value of stretching is overpassed (jDanielsl . Il945l ). The stress carried by fibers 
that break is distributed to the intact fibers of the system, and this gener- 
ates a feedback process that can produce the total failure of the system. The 
description of creep phenomena by the use of fiber bundle model requires 
the consideration of some time dependent process in the model. In some 
cases, it is assumed that the external load increases with time, usually this 
increase is thou ght to be quasi s tatic. This is the "static" FBM considered 
for instance by iPradhan at al. In a slightly different approach in- 

troduced by iKun et al.l (j2003bl ). the external load is kept constant, but the 
system is endowed with some kind of viscoelastic behavior that produce a 
time dependent response. 

An additional effort to set up a model that describes the ful l phenomenol- 



ogy of the three regimes of creep behavior has been made by iNechad at al 
fl2005al lbh. The authors show how previous versions of FBMs fail to obtain 
realistic strain rate dependencies in the Andrade and tertiary creep regime. 
They construct upon the previous models by modifying the rheological prop- 
erties of the fibers, in order to reproduce these regimes. The final conclusion 
is that an Eyring rheology is able to predict a power law in the Andrade and 
tertiary creep regimes. The crossover time t m found with this model is equal 
to tf/2, and the dependence of time to failure with the excess stress a — o c 
is tf ~ (cr — cr c ) -1 / 2 . 

Despite the fact that the modeling presented by INechad at al.l (l2005ai lbl) 
is able to reproduce many of the observed properties of the creep behavior 
of the material, I would like to present here an alternative model that also 
explains many characteristics of the creep process. So first of all it is neces- 
sary to explain why w e need another model of creep rupture. The model of 



Nechad at al.l ( )2005al Jbl) relies on the assumption that fibers respond through 
an ad hoc non-linear rheology. Although the assumed non-linear behavior 
of each fiber can be considered realistic, as it is compatible with macro- 
scopic phenomenology of thermal creep, to apply this macroscopic behavior 
to presumably microscopic elements leaves the impression that we have not 
reached a really fundame ntal description of creep. In ddition, this model (as 
the viscoelastic model of Kun et al. (l2003bh produces a temporal breaking 
sequence of fibers in a rigorous one-by-one process, i.e., avalanche s of ma ny 
breakings do not occur. This forces the authors in INechad at al.l (l2005bU al) 
not to consid er the fitting of acoustic emi s sion data. I n this respect, the model 
presented by iBaxevanis and Katsaounid ( 120071 . 120081 ) seems a bit more realis- 
tic, however it displays some non-physical behavior also (see below). Another 
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additional minor points is that the model in iNechad at al.l (j2005bU al) is not 
able to explain other robust features of the experiments, such as the already 
mentioned constant ratio t m /tf ~ 0.6. 

The second set of reasons to propose a different model for creep rupture 
comes from the fact that the 'new' model I will present, is not actually new, 
nor devised specially to treat this problem, but is a model that has proved 
very successful in reproducing the characteristics of temporal sequences of 
seismic events , giving; in particular a very realistic description of aftershocks 
f|jadal . l2010a] lbh. With rather direct modifications, the modeling of the seis- 
mic phenomena can be adapted to model creep rupture with reasonably good 
results. As we will see, except for the possibility of breaking of elements, the 
model is linear, yet is able to reproduce much of the phenomenology that 
has required the consideration of non-linear rheology in previous attempts. 
This is an important conceptual improvement that justifies the consideration 
of this kind of model. The model is presented in the next Section. Results 
are contained in Section III, and Section IV contains a summary and the 
conclusions. 



2. Model 



The mode l that will be presented here is based on the one that was 
introduced in IJaglal fcOlOaf ). and that has been shown to produce realistic 
sequences of seismic events that quantitatively compare very well with actual 
seismic sequences. Here I will rephrase the main ingredients of the model, 
focusing on its application in the context of creep rupture. I consider a 
set of Nq linear elastic fibers, characterized by the value of their stiffness 
ki (i = 1,...,Nq). The fibers are considered to be equally loaded, i.e., the 
same stretching e is applied to all of them. Fibers break when a maximum 
stretching value ef 1 is reached. Note that for the case of uniform stiffness, 
namely fcj = k Q , the model corresponds to the usual fiber bundle model in 
the global load shearing (GLS) scheme (also known as local load shea r ing;, o r 
democratic fiber bundle) as discussed extensively bv lPradhan at al.l ( 120101 ). 
for instance. Here I will modify the model presented by I Jaglal ( 12010al ) in the 
following way. Let us suppose first of all that the fibers have all the same 
breaking strain e 0; but that they are not anchored to the point of coordinate 
0, but instead to some point of coordinate «j. This means that effectively, the 
stretching that has to be applied to break a particular fiber is e = Eq + Ui. Let 
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us consider now the mechanical energy stored in the system at some given 
value of e. This is given by 

£ = £'|(^) 2 (i) 

i 

where the prime in the sum means that this is extended only to intact fibers, 
namely those with (e — Uj) < sq. The key new point in the system dynamics 
is to consider that the values of ui are not constant but they can evolve 
in time. To work on a safe physical basis, the evolution of u,i will be such 
that the mechanical energy E tends to be minimized in the process. The 
simplest dynamics one can assume for such a temporal variation consists in 
calculating dui/dt as 

^ = -A^ = Xfi = XUe - m) (2) 
dt 5ui J K ' K J 

for intact fibers, where I note the force on fiber % as fa. This would correspond 

to the modeling of each fiber as a linear combination of a harmonic spring 

(with length e — Uj) plus a linear dashpot represented by the variable itj. This 

situation however is not appropriate for the case we are trying to model, since 

in particular, it implies that for any fixed stretching, eventually the system 

reaches a configuration of zero stressQ This would correspond to a viscous 

fluid rather than to the creeping material we are trying to model. Instead of 

the previous situation in which the temporal variation of Ui is proportional 

to the force on the fiber % alone, I propose to model creep assuming that the 

temporal variation of it, will be proportional to the fluctuation of the force in 

the neighborhood of the fiber i. In the case all fibers are intact, it is direct to 

make explicit this prescription in the form of a conserving dynamics, namely 

§=*<r-™=B/.-« (3) 

3 

where the index j labels the four neighbors (on a square lattice) of the site 
i. Some care is however necessary when there are broken fibers, as they have 



1 Note that in the present form this model is similar with the one used by 
Baxevanis and Katsaounia ( 2007 . 20081) . with the difference that they consider that fiber 
rupture occurs at a maximum value of (in the present naming of variables) e + u». In their 
case the unrealistic behavior manifests in the fact that the system eventually breaks for 
any value of applied stress, no matter how small. 
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in some sense come out of the game, and should not contribute to determine 
the temporal evolution of «j. The way in which this is done is the following. 
For each intact fiber i, the temporal variation of Ui is given by 



where the index j is now understood to run only onto intact neighbor fibers 
of the site i. This is the proposed form of the time evolution of Ui that I 
will study in detail in this work. A couple of things are worth to be noted. 
For a given value of e, relaxation through Eq. (j3J) leads eventually to a 
situation in which the force onto any connected patch of intact fibers (if 
any) becomes uniform. This is the relaxed configuration that the system 
approaches, asymptotically. Also, note that if all ki were equal, starting 
from a situation in which all Ui — 0, the previous evolution gives no time 
dependence of U{. This points to the crucial role played by the fact that in 
the present model the values of ki are widely distributed. 

I will concentrate mainly in the case in which an external constant force 
/ per fiber is applied to the system. When one fiber breaks, its load is 
immediately distributed onte the intact fibers, and this can produce the 
immediate failure of some other fibers, generating an avalanche. Whether 
there is a stationary condition in which some of the fibers remain intact, or 
all fibers are eventually broken, depends on the applied load onto the system. 
The protocol to obtain the time evolution of the system, in particular, the 
sequence of fiber breaking is the following. Previously to the application of 
the force, it is assumed that the system had plenty of time to relax, in such 
a way that the forces onto each fiber have become equal to zero everywhere. 
This means that all ui = 0. At t = a nominal force / per fiber is applied. 
Assuming that initially all fibers are intact, the following sequence of steps 
is performed (primed sums run only onto intact fibers): 

1- The value of e is calculated from Yl'i( £ ~ u i)h = fN 

2- Check if any fiber has overpassed the breaking threshold, i.e., if (e — 
Ui) > Eq. If yes, consider this/these fiber/s as broken, and return to 1. 

3- Once no additional broken fibers occur in step 2, increase time in a 
small time step dt, calculate the new U{ values using (BJ, and return to step 
1. 

In this way, given an externally imposed value of the force /, the time 
dependence of the stretching e(t) can be followed. e(t) will typically have 




(4) 
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periods of continuous evolution when no fibers breaks, and abrupt jumps 
occurring when fibers break. I emphasize that as the force carried by any 
fiber that breaks is distributed onto the intact fibers (through steps 1 and 
2), a cascade of events can occur in which a finite number of fibers break 
exactly at the same time. This i s an interesting variant with respect to o ther 
viscoelastic fiber bundle model (INechad at all l2005bt iKun et all l2003bl ). in 
which the breaking of fibers occurs in a strictly one-by-one sequence. 

A technical comment on the implementation of the previous rules is ap- 
propriate here. If we want to follow in full detail the bursts of breaking 
activity in the system, the time step dt in step 3 above has to be tuned in 
such a way that it never produces by itself the breaking of more than one 
fiber at a time. This is easily done as the precise value dto necessary to break 
the next fiber can be calculated due to the linearity of the springs. Then in 
this approach, the value of dt is chosen dynamically at each passage through 
step 3. This full detail process however, is very time consuming if we are 
only interested in the overall evolution of e with time. In this fixed 
value of dt is used, and in this time interval a few different fibers can break 
independently, however, the overall evolution of e(t) is very well reproduced. 
Now I will present detailed results on the numerical implementation of this 
model. 



3. Results 

3.1. Results for average quantities 

I already indicated that a broad distribution of stiffness ki is necessary 
for the model to be sensible. The precise form of this distribution is however 
not crucial for the observed properties. I will present results for a uniform 
distribution of k between and k max = 1. Additional calculations using a 
Weilbull distribution have shown no appreciable qualitative differences. 

Depending on the value of the external load per fiber /, the configura- 
tion of the system at a very large time corresponds to one of two different 
possibilities. If / is lower than some critical value f c only a fraction of fibers 
break, and the value of e remains bounded as t — > oo, whereas if / > f c , all 
fibers break in a finite time. 

In Fig. Ufa) we see the evolution of the stretching e as a function of time, 
for a single realization in a system of 1000 x 1000 fibers. The qualitative 
difference between the two regimes f > f c and / < f c , is clearly observable. 
For the case / > f c , the time to failure in the system diverges as (/ — f c )~ a , 
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Figure 1: (a) Strain as a function of time at different values of the applied load in a single 
realization on a system of size 1000 x 1000. In the inset, the time to failure tf for f > f c 
is plotted as a function of load, showing a power law divergence approaching f c — 0.2763. 
Points are the numerical results, dotted line has a slope of -1.25. (b) The approach of e 
to its asymptotic value for the case / < f c is shown, and is seen to follow a power law 
behavior (dotted lines have slopes -0.45 and -0.55, for reference). The best determined 
asymptotic values of the strain as a function of the applied force are plotted in the inset, 
where the critical force f c is marked as a vertical dotted line. 



with a ~ 1.25 ±0.01, as can be seen in the inset. For the present realization, 
f c turned out to be f c = 0.2763. This value shows rather large fluctuations in 
different realizations. In panel (b) of Fig. (TJ I show the approaching of e to the 
asymptotic values, where a power law dependence is clear. The asymptotic 
values have been adjusted to obtain the best power law fitting, and they are 
shown as a function of / in the inset. It is seen that they follow a smooth 
behavior as f c is approached, without obvious signs of any singular behavior. 
The power law approach of e to the asymptotic values is compatible with 
the Andrade creep law. Note however that the exponent ~ 0.5 corresponds 
to a dependence of the strain rate e of the form e ~ 1/£ L5 , whereas the 
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Figure 2: Strain rates e = de/dt as a function of time normalized to time-to- failure 
in the run away regime (/ > f c ). The data shown correspond to an average over 120 
realizations in systems of size 1000 x 1000. In first panel time scale is linear, and the time 
of minimum strain rate has been emphasized. In the second panel time scale is logarithmic, 
to display the power law decay in the primary creep regime (note however the larger-than- 
one exponent compared to usual Andrade creep). In the last panel a logarithmic time scale 
approaching tf is used, to emphasize the power law increase of strain rate as approaching 
failure. 
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usual exponent for strain rate in the Andrade regime is about 2/3. Note 
however that the usual Andrade regime for / — f c has eventually to cross 
over to a different dependence that guarantees a finite strain at infinite time, 
(otherwise the lower-than-one exponent would produce a divergent strain). 
Here the larger-than-one exponent obtained gives a finite strain at infinite 
time. 



To m ake a better comparison with the results presented by lNechad at al 



( 12005bl ). in Fig. [2] I present results for / > f c , showing the values of strain 
rate in the system e = de/dt, as a function of t/tf, under different scaling of 
the temporal axis. In the first panel the temporal axis is linear, and rescaled 
with the time to failure tf. Here we see that e displays qualitatively the three 
mayor creeps stage before rupture: a primary regime of decreasing creep rate, 
a second regime of almost constant creep, and a tertiary regime of accelerated 
creep towards failure. 

From this figure we see that the time at which the creep rate is minimum 
is between at 0.6-0.7 of tf, independently of /, a result that was observed 
to occur systematically in real samples, and that was not captured by the 
models available up to now. The logarithmic temporal scale in panel (b) of 
Fig. [2j shows that the first stage of creep in the model corresponds indeed to 
an Andrade regime of potential decay of creep rate with time. The exponent 
observed for this decay is seen to be about 1.3, which is larger than the 
usually observed exponent that is typically near 2/3. The logarithmic plot 
with respect to (tf — t) in panel (c) contains a very interesting outcome of 
the model: strain rate shows a tendency to diverge at tf as (tf — i) -1 , which 
is in fact the result that is observed in a variety of experimental cases. It has 
to be reminded here, that fiber bundle models provi de typically a divergen ce 
of strain rate as (tf — t)^ 1 ^ 2 , whereas the model by iNechad at al.l (2005bJJaj) 
obtains a systematic (tf — t)^ 1 dependence relying on a specific non-linear 
choice for the rheology of the fibers. 

An interesting additional point on the "divergence" of e at tf has to be 
noticed. Although the plot in Fig. [2](c) shows an increase of e by about two 
orders of magnitude compatible with a (tf — t)^ 1 dependence, this behavior 
cannot be extrapolated to t — > tf. In fact, if this was the case, then the 
value of e itself would diverge at tf but this is nonsense: there are no fibers 
that can withstand an arbitrarily large deformation. The situation is that e 
increases as (tf —t)^ 1 (and eas 1/ log(i/ — t)) up to some time slightly before 
tf, at which a macroscopic event breaks all remaining fibers (which is a finite 
fraction of the total system) all of the sudden. We will confirm this scenario 
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in the next section where I concentrate on the study of bursts of activity in 
the system. 

3.2. Burst sizes distribution 

The results in Figs. [1] and [2] focus on the average behavior of the system, 
namely, the evolution of the stretching e or its time derivative e on time, 
and are indications that the model is able to capture in a quite realistic 
way the time evolution of strain that is observed in materials experiencing 
creep rupture. Now we will see that fluctuation in these quantities, asso- 
ciated in particular to the existence of discrete breaking events, reveal also 
very interesting behaviors that can be compared for instance with results of 
experiments on acoustic emission. 




Figure 3: Main upper panel: A sequence of all breaking events in a system of size 1000 
x 1000, for a load slightly larger than critical (f / f c — 1-05). The average burst size A 
as a function of time is indicated in the inset, in log scale (the power law of the straight 
line has an exponent 0.25). The cumulative number of events E(t), and the cumulative 
number of broken hbers N(t) as a function of time is shown in the panel below. 



Some previous work have focused on the size distribution of bursts. How- 
ever, in some other models the burst do not occur exactly as a function 
of time under a constant force as in our case. In the much studied static 
model considered in particular by IPradhan at al.l (120 lOl ). it is assumed that 



the external force increases quasi-statically in time, then the bursts occur 
because / is increasing, and not directly because of the existence of creep 
phenomena. Anyway this is a case that has been studied in much detail, and 
is the one with which we should compare, also because there are interest- 
ing analytical predictions for the size distribution of bursts in this case. In 
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fact, it is known that if bursts are counted during the whole load process, 
from / = to complete failure, the number of bursts of size A follows a 
power law of the form A ~ A - ^, with £ = 5/2. If however, the size of the 
bursts is collected only in a small time interval before breaking occurs, the 
ex ponent predicted is £ = 3/2. Some other viscoelastic models, as the one 
by iKun et al. fl2003bl ) produce a one-by-one breaking of fibers, and in that 
case a non-trivial size distr i bution of events cannot be defined. The model of 
Baxevanis and Katsaounisl (l2007h seems to be more realistic in this respect, 



and the £ = 5/2 and £ = 3/2 power laws for the total and close-to- rupture 
avalanches have been obtained in this case also. 

The model we are studying here gives a non trivial time sequence of 
bursts for any fixed vale of /. In Fig. [3] we see a typical time sequence of 
all breaking events in the system, for a value of / slightly above the critical 
value f c . The cumulative number of broken fibers N(t) is qualitatively very 
similar to the evolution of e(t) discussed in the previous section. This number 
is seen to approach a fraction of about 0.55 of N as t approaches tf. At this 
point, a macroscopic breaking event occurs, in which the rest of the system 
breaks. Note that the fraction of fibers that breaks in this single event (about 
0.45iVo) is rather independent of the value of N itself, and on the value of / 
(assumed larger than f c ). The existence of this macroscopic breaking event 
is reasonable in the present model, where the breaking of a fiber implies the 
immediate load increase on all the other fibers. In a more realistic situation, 
particularly as the time of global failure is approaching, even very small 
delays in the process of redistribution of stresses will have an important 
influence on the dynamics of failure. I expect that this consideration results 
in a continuous evolution of N(t) towards Nq at tf in that case. The results 
for the the average size A of the bursts as a function of time (inset in Fig. EJ) 
shows that A increases as a power law as tf is approached, namely A(t) ~ 
\/{t-t f ) z , with z~ 0.25. 

The distribution of bursts corresponding to the sequence in Fig. [3] is 
shown in Fig. H] for different time interval [to>^/] with < t < tf. The 
complete burst distribution (to = 0) shows a consistent power law with a 
decaying exponent around 2.8. As to increases, a crossover and a coexis- 
tence of two different power laws is observed. For to — > tf the distribution 
tends to be described by a single power law with an exponent close to 1.5. 
These results are qualitatively similar (except for a sl ightly larger global ex- 
pon ent) to those obtained in the mo dels studied by Pradhan at al. f 2010l ) 
and iBaxevanis and Katsaounisl (120081 ). 
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Figure 4: Distribution of bursts for the sequence in figure [3l Results are presented for dif- 
ferent time intervals, namely [to>i/] for different values of to (tf — 80.82 for this particular 
realization). The distribution is seen to pass from a power law with a decaying exponent 
~ 2.8 for the full distribution (i.e., to — 0), to a 1.5 power law when to —¥tj. 



Results for bursts sizes can be cast is a slightly different form in order to 
compare with available results on acoustic emission. In particular, one of the 
magnitudes that has been experimentally recorded is the cumulative energy 
released as a function of t ime, on approaching failure. It has been obtained 



(jGarcimartin at al.l . 119971 ; iGuarino et all 120021 ) that this magnitude shows a 



power law increase on approaching failure with a dependence of the form ~ 
(tf — t) 7 , with reported values for 7 close to 0.25. The experimental detection 
and processing of an acoustic signal involves some non-trivial steps that make 
the link with other physical properties of the system not a straightforward 
procedure. However, for the present model we can reasonably consider that 
a burst involving the breaking of A fibers produces an acoustic emission 
signal of energy ~ A 2 , as the energy is experiment ally defined as the time 



integral of th e squar e of the amplitude detected (jYagamuchi at al.l . 11991 



Nechad at al.l . l2005al ). Also, we should consider that events involving less 



that some minimum value A of broken fibers are not detected because lack 
of sensitivity. In Fig. [5] I show the result for the accumulated energy released, 
calculated by integrating in time the values of A 2 , for events with A larger 
than a threshold Ao, corresponding to the sequence of events shown in Fig. 
[3j The results show the right tendency to diverge as t — > tf. The definition of 
a well defined power law ~ (tf — /)~ 7 is however strongly dependent on the 
threshold A that is used. Values of 7 between 0.2 and 0.5 can be obtained 
by changing A . 
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Figure 5: Cumulative energy released as would be detected in an acoustic emission experi- 
ment corresponding to the sequence of events shown in Fig. [3J Different curves correspond 
to different detectability threshold A , as indicated. A tendency to increase as a power 
law as tf is approached is observed, however, the exponent is dependent on the threshold 
used. Two lines with exponents 0.2 and 0.5 are drawn for comparison. 

3.3. Interevent times 

Another quantity that is expected to provide information on microscopic 
processes in the system, and that is experimentally accessible in acous- 
tic emission measurements is the distribution of inter event times (IETs), 
namely, the time intervals St between successive breaking events in the sys- 
tem. For systems in which bursts occur in an uncorrelated time sequence, an 
exponential (Poissonian) distribution of IETs is expected. Deviations from 
the exponential distribution might be taken as signs of non-trivial correla- 
tions in the activity of the system. However, this interpretation needs to 
be used with care, at least. In fact, in Fig. [6ta) we see the full sequence 
of and IETs for of an applied load slightly below the critical value. 

The activity in the system decays in time towards a stationary state, and 
this is seen in the systematic increase of IET with time. The distribution if 
IETs for this case is shown in part (b) of Fig. |HJ It is seen there that IETs 
show a power law distribution over more than seven orders of magnitude. 
The origin of this power law is however trivially related to the fact that the 
average IET as a function of time is itself a power law increasing function 
of time. By calculating the distributions of IETs for partial time intervals 
(shown in Fig. El^b)) makes clear that the power law appears only as a result 
of the integration of different time windows. 

It is thus evident that if non-trivial correlations between events exist, 
their signatures have to be searched within small time windows, in which 
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Figure 6: (a) A sequence of interevent times St in a single realization in a 1000 x 1000 
system, for a load slightly below the critical value. The average interevent time increases 
approximately as t 1 - 5 (dotted line), (b) Full symbols: distribution of interevent times 
in the complete time interval [0,oo]. Open symbols: distribution of interevent times for 
different time intervals ti-t%, as indicated. 



the process can be considered stationary. This fact seems not to have been 
taken into account properly in some previous studies of IETs distributions. 
In Fig. [7] I show results in a system of 400 x 400 for a force slightly above the 
critical value. A smaller system size is used in this it makes easier to 

detect correlations between events. In order to have a reasonable statistics, I 
run different realizations, and the results have been combined by scaling the 
occurrence time of the events by the corresponding failure time tf of that 
particular realization. The results for the distribution of IETs calculated for 
small time windows (Fig. [7(a)), show a clear departure from the Poissonian 
behavior. The curvature of the data indicates an increased probability of 
events to occur at small or large time intervals, whereas probability for inter- 
mediate time intervals is depressed. This indicates a tendency of the data to 
cluster in time, which is also pretty clear in the individual realization shown 
in the lower part of the figure. Temporal clustering of events in the present 
model is reminiscent of temporal clustering of earthquakes, in the form of 
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aftershocks. In the present model, temporal clustering of events is a direct 
consequence of the relaxation mechanism present in the model, as it was 
qualitatively expla i ned in the version of the model devised to study seismic 
phenomena (Jagla, 2010al bl). This kind of clustering should be detectable in 
acoustic emission experiment. 




Figure 7: (a) Distribution of IETs within small time windows [t\/tf, £2/*/] a s indicated, in 
systems of 400 x 400, for / slightly abobe f c (the result of an average over 120 realizations 
is shown). The data do not follow an exponential distribution, as systematic deviations 
of a straight line plot are observed. This indicates temporal clustering of the events, (a) 
A particular sequence of events and IETs where the temporal clustering is observed. A 
temporal region of enhanced (depressed) activity is highlighted by the continuous (dotted) 
arrow. 



4. Conclusions 

In the present paper a model to describe the creep rupture of materials 
has been presented. The model is in the class of viscoelastic fiber bundle 
models, and the viscoelastic relaxation mechanism considered involves a ten- 
dency to locally uniformize the stresses within the system. The model is a 
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slight modification of one that has been proved useful to explain many char- 
acteristics of the sequences of events (earthquakes) in a seismic context. The 
only non-linear ingredient of the model is the possibility that fibers break 
when a maximum stretching is overpassed, in particular, I do not rely upon 
any non-linear rheology of the elements of the system. Using this model, I 
have obtained a description of the creep rupture process that in a few re- 
spects is an improvement with respect to the results in previous modeling. 
In particular, a regime of primary (Andrade) creep has been obtained, in 
which the strain rate has a power law dependence. The exponent of this 
dependence, is however larger than one, contrary to the usually observed 
exponent around 2/3. A region of accelerated creep preceding rupture was 
obtain, in which the strain rate increases as ~ (t f — t) -1 , where tf is the 
time of failure. The time of minimum strain rate is systematically observed 
to occur at about 0.6 of the time of failure, coincident with experiments. 
Bursts sizes distributions around rupture with an exponent close to 3/2, and 
for the whole temporal sequence with an exponent slightly above 5/2 where 
obtained, much similar to results with previous modeling. The accumulated 
acoustic emission energy in the model was shown to increase as a power law 
as the time of failure is approached. The observed exponent is however de- 
pendent of the cut off used as a threshold for detecting events, and varies 
between 0.2 and 0.5. A power law distribution of interevent times (IETs) is 
trivially originated in a temporal variation of the average IET. Analysis of 
IETs in time windows of constant average activity reveal that the activity is 
correlated to some extent, displaying temporal clustering of the events. This 
was indicated to be reminiscent of the aftershocks that are produced with 
this model when used in the context of seismic phenomena. 
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